Magnetic Braking and Protostellar Disk Formation: 

Ambipolar Diffusion 

Richard R. Mellon 1 & Zhi-Yun Lr 1 
ABSTRACT 

It is established that the formation of rotationally supported disks during 
the main accretion phase of star formation is suppressed by a moderately strong 
magnetic field in the ideal MHD limit. Non-ideal MHD effects are expected 
to weaken the magnetic braking, perhaps allowing the disk to reappear. We 
concentrate on one such effect, ambipolar diffusion, which enables the field lines 
to slip relative to the bulk neutral matter. We find that the slippage does not 
sufficiently weaken the braking to allow rotationally supported disks to form for 
realistic levels of cloud magnetization and cosmic ray ionization rate; in some 
cases, the magnetic braking is even enhanced. Only in dense cores with both 
exceptionally weak fields and unreasonably low ionization rate do such disks start 
to form in our simulations. We conclude that additional processes, such as Ohmic 
dissipation or Hall effect, are needed to enable disk formation. Alternatively, the 
disk may form at late times when the massive envelope that anchors the magnetic 
brake is dissipated, perhaps by a protostellar wind. 

Subject headings: accretion disks — ISM: molecular clouds and magnetic fields 
— MHD — stars: formation 



Introduction 



Disk formation is an integral part of star formation that has been studied f or a long time. 



Early works concentrated o n the collapse of rotating, non-magnetic cores (e.g-. lTerebey et al. 



19841 ) . iBodenheimerl (1995T ) reviewed these works, and noted a number of unresolved prob- 
lems (see also lBosslll998l ). Topping the list is the effect of magnetic braking on disk formation. 



The magnetic properties of star-forming dense cores are reasonably well constrained. 
There is now ample evidence for ordered magnetic fields on the core scale from polarization 
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of dust emission (e.g., I Ward- Thompson et al.ll2000l ). A spectacular recent example is the 



AU scale ( 


Girart et a 


2006) 


Troland & Crutcher 


( 


2008h < 



2006). The magnitude of the magnetic 



survey of dark cloud cores in OH, an d determined a mean mass- to-flux ratio of A ~ 4.8 in 



units of the critical value (2nG 1 ^ 2 )^ 1 (jNakano fc Nakamuralll978l ); correction for uncertain 
projection effects may bring this value to ~ 2. The implication is that dense cores are 
typically moderately strongly magnetized, with a dimensionless mass-to-flux ratio of a few 
to several. A relatively low value of A ~ 2 — 3 is expected for dense cores formed out of 
magnetically subcritical clouds (with A < 1) through ambipolar diffusion (e.g . , iLizano &: Shu 
19891 : iBasu fc Mouschoviasl ^994 : iNakamura fellbood : Adams fc Sh~u1l2007h . The values of 
A are expected to have a wider spread for cores formed out of turbulent compression. Nev- 
ertheless, the cores tend to be more magnetized relative to their masses than the cloud as 
a whole, because only a fract ion of the mass along a flux tube tha t threads the cloud is 
compressed into the core (e.g.. lDib et al.ll2007l ; iTilley fc Pudrita 120051 ). It is unlikely for the 
cores to have a value of A more than several, unless the cloud as a whole is magnetized to 
an unrealistically low level. 

There have been a number of recent calculations that included magnetic fields in the 
collapse of rotating cores, focusing mostly on the early phase before the mass of the central ob- 
ject reaches the stellar ra nge (e.g.. lTomisakalll998l : iMachida et al.ll2006l : iBanerjee fc Pudritz 
20061 ; iPrice &: Bate! 120071 ). Other studies have concentra ted on the angular momentum evo- 
lution and disk fo rmation in the mai n accretion phase ([Allen et al~ 2003 ; Galli et al. 2006 



Mellon fc Li 2008 . hereafter Paper I; Hennebelle fc Fromang! 20081 ). These studies find that 
magnetic braking efficiently removes angular momentum from the collapsing matter, pre- 
venting a large, 10 2 -AU scale rotationally supported disk from forming in realistic cores in 
the ideal MHD limit. However, such disks are routinely observed around young stars, at 
least in relatively late, Class I and II phases. How can the disk be saved? 

Molecular cloud cores are lightly ionized, so perfect coupling between the magnetic field 
and matter assumed in ideal MHD is not expected. As the density increases, the coupling 
is weakened first by ambipolar diffusion, then the Hall Effect, then Ohmic dissipation. In 
this paper we will concentrate on ambipolar diffusi on as a first step. The effec t of ambipolar 
diffusion on magnetic braking in the formation (IBasu &; Mouschoviasl 1 19941 ) and collapse 
(IKrasnopolsky fc Konigll 120021) of rotating cores has b een examined previously using the 
thin disk approximation. IKrasnopolsky fc Konigll (120021 ) parametrized the braking strength 
and found that disk formation can be suppressed in the ambipolar diffusion limit when the 
braking parameter is large. We use two-dimensional models which allow direct calculation 
of the strength of the magnetic braking. We find that ambipolar diffusion does not weaken 
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the magnetic braking enough to allow rotationally supported disks to form under realistic 
conditions. The implication is that additional processes must be found to save the disk; two 
possibilities are discussed in section 5. 



Model Formulation 



We study disk formation and evolution during the main accretion phase after a central 
object has formed. The problem setup is identical to that of Paper I (to which we r e fer th e 
reader for details), except for the inclusion of ambipolar diffusion. As in I Allen et al.l (120031 ). 
we idealize the initial configuration of the main accretion phase as a rotating self-similar 
singular isothermal toroid supported against self-gravity partially by thermal pressure and 
partially by magnetic fields (see figure 1 of Paper I for an example). The toroid is char- 
acterized by three parameters, the isothermal sound speed (c a ), mass-to-flux ratio (A) and 
rotational speed (v ), which is taken to be spatially constant to preserve the self-similarity 
of the configuration. 

Our treatment of ambipolar diffusion follows that of Shu (1991), where the magnetic 
field is assumed to be tied to ions, which drift relative to neutrals with a velocity 



1 



(V x B) x B 



T. 



where 7 = 3.5 x 10 13 cm 3 g _1 s" 1 is the ion-neutral drag coefficient and p n (pi) and v n (vi) 
are the density and velocity of the neutrals (ions), respectively. We adopt an ion density of 
the form 

P, = C pT 

with the coefficient 



(2) 



C = 3 x 10~ 16 



C 



10- 17 s- 



1/2 



cm ~ 3 / 2 g 1 ' 2 



(3) 



where ( is the cosmic ray ionization rate (lElmegreenl Il979l ). This simple form has the 
advantage of preserving the self-similarity of the collapse, which provides a powerful check on 
our numeric ally obtained soluti ons. We implement the ambipolar di ffusion into the ZeusMP 
MHD code (IHayes et al.l 120061 ) using the fully explicit method of iMac Low et al.l (119951 ). 
As is well known, for such explicit methods to be stable, the limiting timestep must be 
proportional to the grid size squared, which puts stringent constraints on spatial resolution. 

We solve the governing MHD equations including ambipolar diffusion in a spherical 
coordinate system (r, 6, <ft) assuming axisymmetry. The grid spacing is logarithmic in radius 
r and constant in polar angle 9, with a computational domain extending radially from 10 14 



-4- 



to 2 x 10 17 cm and to rr in angle. There are 120 grid points in the radial direction, and 
60 angular points. The smallest grid size is ~ 10 13 cm, near the inner boundary. The 
standard outflow condition is imposed at the outer radial boundary. The inner boundary 
is a modified outflow boundary with the mass accreted across the boundary added to the 
central point mass and a torque free condition imposed on the magnetic field (i.e., = 0, 
see Paper I). We adopt the oft-used broken-power law equation of state that is isothermal 
below p = 10~ 13 g cm~ 3 and adiabatic with 7 = 7/5 above. 



3. Standard Model 



We first illustrate the effects of ambipolar diffusion on magnetic braking and disk forma- 
tion using a model with a fiducial cosmic ray ionization rate ( = 10~ 17 s -1 and a mass-to-f lux 
ratio A = 4; the latter is close to the mean value inferred by iTroland fc Crutcherl (120081 ) in 



their Arecibo OH Zeeman observations of dark cloud cores. Other combinations of ( and A 
are considered in the next section. For all models, we set the sound speed to c s = 0.3 km/s 
and the rotational speed to v = c s /2 = 0.15 km/s, which corr esponds to an angular speed 
of 3 km s _1 pc _1 , typical of dense cores (IGoodman et al.lll993l ). 



A quantity of primary importance to star formation is the mass accretion rate, M. In 
Fig. [TJ we show M as a function of time. In our standard model, M converges to a constant 
value, after a short period of adjustment. The initial deviation from the constant is due to 
the zero infall velocity at t — and the small point mass that we put at the center to induce 
the collapse. The accretion rate quickly approaches a constant value indicating the solution 
has reach ed the exp ected self similar state at late times. The converged value of M is close 
to that of lShul (119771 ) for a singular isothermal sphere. One might naively expect the rotation 
to retard the collapse significantly, leading to a reduced mass accretion rate onto the central 
object. Apparently, the magnetic braking is strong enough that rotation is not a significant 
barrier to mass accretion, even in the presence of ambipolar diffusion. The accretion rate is 
comparable to that in the ideal MHD case (Paper I), suggesting that the strength of braking 
is not significantly reduced by ambipolar diffusion. The accretion is highly episodic in the 
ideal MHD case, but is steady in the presence of ambipolar diffusion. 

Fig. [2] shows a snapshot of the standard model at a representative time of 5.85 x 10 11 sec. 
As in the singular isothermal sphere, the collapse occurs inside out, with the bulk of the 
envelope material beyond the radius of c s t = 2 x 10 16 cm remaining nearly static. Inside this 
radius, the matter accelerates towards the center. The col lapsing flow is defl ected towards 
the equator by pinched field lines, forming a pseudodisk (jGalli &: Shulll993l ). As the flow 



collapses, it is expected to spin up due to conservation of angular momentum. However, the 
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rotation speed tends to decrease as the collapsing material moves inward, especially in the 
equatorial region, indicating significant loss of angular momentum due to magnetic braking. 

To illustrate the braking more quantitatively, we plot the equatorial velocity profiles 
in Fig. [31 As one moves from large radii to the center we observe four distinct regions 
in the radial velocity profile (top panel). Outside of r = 2 x 10 16 cm the matter remains 
nearly static. In the region 6 x 10 15 cm<r<2 x 10 16 cm the matter accelerates toward 
the center as gravity dominates the magnetic forces, since the mass-to-flux ratio of the 
collapsing core is significantly greater than unity. The collapsing material is decelerated 
in the region 3 x 10 15 cm<r<6 x 10 15 cm, where the outward magnetic forces exceed the 
inward gravitational pull. This region does not exist in the ideal MHD limit; its presence 
is due entirely to ambipolar diffusion, which enables the magnetic field lines that would 
have be en dragged into the central object in the ideal MHD limit to pile up outside the 
object Li fc McKeel Il996l : Iciolek fc Korng] Il998h . In this decelerating matter, the ion- 



neutral drift velocity becomes large because of a large Lorentz force. At r ~ 3 x 10 15 cm 
the neutral matter begins accelerating towards the center again. However, the ions remain 
nearly static at r ~ 1.5 x 10 15 cm, indicating that the magnetic field lines barely move in that 
region. The deceleration region has th e same physical origin as the hydromagnet ic accretion 



shock induced by ambipolar diffusion (jLi fc McKeel Il996l ; ICiolek &: Konigllll998l ). However 



since the gas outside of the decelerating region is collapsing subsonically, a hydromagnetic 
accretion front rather than a shock is formed. At r<10 15 cm the ions (and magnetic field) 
begin reaccelerating towards the center. The rapid collapse of both ions and neutrals at the 
smallest distances from the star clearly suggest the absence of a rotationally supported disk. 

From the lower panel of Fig. [3j it is clear that the rotation speed v</, does not increase 
with decreasing radius, which is opposite of what is expected based on angular momentum 
conservation. For example, in the region of accelerating collapse (6 x 10 15 cm<r<2 x 10 16 cm) , 
remains constant while the radius decreases by a factor of ~ 3. This indicates significant 
loss of angular momentum. This loss can be understood by examining the ratio of toroidal 
to poloidal magnetic fields. We find that the magnetic field lines get more twisted as the 
collapse accelerates inwards. When deceleration occurs, compression increases the poloidal 
and toroidal field strength, which leads to a stronger braking that is responsible for the steep 
drop in rotational speed in the deceleration region (the shaded region in Fig. [3]). By the time 
the fluid reaccelerates towards the center, there is little angular momentum left in the gas to 
maintain a significant twist in the magnetic field, resulting in near radial infall to the center. 
The supersonic infall velocity and small rotation speed at small radii leaves no doubt that 
rotationally supported disk formation is suppressed in our standard model. 
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4. Effects of Ionization Rate and Magnetic Field Strength 

Our collapse solution depends on the cosmic ray ionization rate, which is somewhat 
uncer tain. For dense cores , the typical values lie in a range between (1 — 5) x 10~ 17 s" 1 
(e.g., iMaret k, Berginl 120071 ). although a wider spread is also possible ( ICaselli et al.lll998l ). 
To examine the sensitivity of our results to the ionization rate, we vary Q by a factor of 10 
in either direction. 



The higher ionization rate case has a larger ionization fraction by a factor of v 10, which 
leads to stronger coupling between the matter and magnetic field. Its collapse is qualitatively 
similar to the standard model, with a few quantitative differences. The stronger coupling 
allows the magnetic fields to be brought closer to the central object, leaving a larger magnetic 
flux in a smaller region. As a result, the deceleration region is closer to the central object 
(see dashed lines of the top panel of Fig. H|). The stronger field and stronger coupling keeps 
the ion radial speed low for a more extended region during the reacceleration. The stronger 
field also results in enhanced braking which removes essentially all the angular momentum at 
small radii (see dashed lines of the second panel, noting that the rotation speeds for ions and 
neutrals are practically indistinguishable). Increasing the ionization rate further suppresses 
rotationally supported disk formation. 



The lower ionization rate yields an ionization fraction smaller by a factor of v 10, which 
leads to weaker coupling between the matter and magnetic field. The ions start to decel- 
erate earlier, due to the pile up of magnetic flux at larger radii. However, the neutrals slip 
through the ions without being appreciably deccelerated due to weak coupling, resulting 
in faster diffusion. This results in a faster collapse of the neutrals and a slower collapse 
of the ions compared to the standard model (see the dot-dashed lines in top panel). The 
increased diffusion requires a smaller ambipolar diffusion timestep, which increases compu- 
tational time3- The weak coupling and reduced magnetic field strength decrease the strength 
of magnetic braking, allowing the neutrals to spin up as they collapse. However, there is still 
enough braking to prevent a rotationally supported disk from forming inside our computa- 
tional domain, as evidenced by the rapid radial collapse at small radii. We conclude that 
for moderately strongly magnetized core of A = 4, disk formation is suppressed for realistic 
values of the cosmic ray ionization rate. 

The mass-to-flux ratio of dense cores is uncertain. Our standard model is based on the 



mean mass-to-flux ratio inferred by iTroland fc Crutcherl (120081 ) , but many cores have only 



1 Although the two models of C = 10~ 18 s _1 in Fig. |4] were evolved for a shorter time (2.95 x 10 11 s 
compared to 5.85 x 10 11 s for other, better coupled models), they still converged to a self similar solution. 
Their velocity profiles have been self similarly scaled for comparison. 
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lower limits on the mass-to-flux ratio. We investigate the collapse model in more weakly 
magnetized cores with A = 13.3 (same as the standard ideal MHD case of Paper I), which 
corresponds to a field strength of 7.35 \xG on the scale of the typical core radius of 0.05 pc. 
We take this case as a lower limit of realistic magnetic field stren gths based on the media n 
field strength of cold neutral HI structures of ~ 6 /iG inferred by lHeiles fc Trolandl (120051 ). 

We again consider three values for the ionization rate ( = 10~ 16 , 10~ 17 , and 10~ 18 s -1 , 
concentrating first on the case with the fiducial value of 10 -17 s _1 . In contrast to the standard 
model with the same £, there are two equatorial decceleration regions in this weaker field 
case instead of one (see solid lines in the third panel of Fig. H]). The accelerating infall 
is first decelerated near r ~ 10 16 cm. This deceleration is not a new feature caused by 
ambipolar diffusion; it was already present in our previous ideal MHD calculation: it is a 
"magnetic barrier" caused by the bunching of magnetic field lines at the interface between 
the collapsing envelope and an expanding, magnetic braking-driven bubble (see figure 3 
of Paper I). Inside the barrier, the material recollapses inwards, at an increasingly high 
speed, until a radius of ~ 10 15 cm, where it is decelerated for a second time. The second 
deceleration is caused by ambipolar diffusion, which allows magnetic flux to pile up at small 
radii, instead of being dragged into the center. As in the standard model, the magnetic 
pressure in the deceleration region becomes comparable to the ram pressure of the collapsing 
material. However, in this model the deceleration region occurs at a smaller radius, allowing 
the infalling gas to collapse supersonically, causing a hydromagnetic shock instead of a front. 
The slowdown of infall compresses the collapsing material, leading to a stronger poloidal 
and toroidal magnetic field that enhance the rate of magnetic braking. The enhancement 
is responsible for the dip on the rotational speed profile (see solid lines in bottom panel). 
Interior to the deceleration region, the collapsing flow accelerates for a third time, spinning 
up as it collapses. Nevertheless, there is insufficient angular momentum left in the material 
for a rotationally supported disk to form before the inner boundary is reached. 

The dynamics of the more ionized, ( = 10~ 16 s _1 case is similar to that of the £ = 
10~ 17 s _1 case down to the magnetic diffusion- induced second deceleration near ~ 10 15 cm. 
Near this region and interior to it, the infall speed is also similar, but the rotational speed 
is drastically different (see dashed lines in the bottom two panels). The field compression in 
the deceleration region increases the rate of magnetic braking to such a degree as to cause 
a counter-rotation in the better coupled case. The better field-matter coupling and stronger 
magnetic field lead to efficient removal of angular momentum at small radii, suppressing disk 
formation completely. 

Only in the most weakly ionized case of £ = 10~ 18 s _1 does a rotationally support 
disk form. The disk is most obvious in the radial velocity curve (dot-dashed lines in the 
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third panel). As in the more strongly coupled cases, there is accelerating collapse at large 
radii. The collapse is only weakly slowed by the magnetic barrier, because the neutrals are 
weakly coupled to the field. The ions and neutrals are reaccelerated, then decelerated for a 
second time. Unlike the more strongly coupled cases, this second deceleration is not due to 
the ambipolar diffusion-enabled pile up of magnetic flux, but is caused by centrifugal force 
from rapid rotation (see dot-dashed lines in the bottom panel). The combination of a weak 
magnetic field and weak matter-field coupling renders the braking too inefficient to remove 
enough angular momentum to suppress disk formation in this case of extreme parameters. 



5. Discussion and Conclusion 



We have studied the collapse of rotating molecular cloud cores magnetized to mod- 
erate degrees, concentrating on a dimensionless mass-t o-flux ratio A = 4, as suggested by 
recent Zeeman observations (ITroland fc Crutcherl 120081 ). In the ideal MHD limit, we have 
previously shown that a slight twist of the field lines in such cores is sufficient to remove 
essentially all of the angular momentum of the collapsing matter and suppress the formation 
of rotat ionally supported disks during the protostellar accretion phase completely (Paper I ; 



see also Galli et al. 



20061 ; iHennebelle &: Fromangll2008l ) . One may expect ambipolar diffusion 



to reduce the efficiency of magnetic braking and potentially enable the disks to form, since 
the twisting of field lines by the rotation of neutrals will be reduced by diffusion. However, 
for realistic cosmic ray ionization rates, the small amount of field twisting needed for an- 
gular momentum removal drives a drift speed between the ions (and the field lines tied to 
them) and neutrals in the azimuthal direction that is much smaller than the neutral rotation 
speed (see the lower panel of Fig. [3]). By itself, the azimuthal slippage of field lines relative 
to neutrals does not reduce the strength of the toroidal magnetic field (and the braking 
rate) significantly. The ambipolar diffusion is expected to modify the poloidal magnetic field 
as well, especially in regions where the magnetic force is large. The modification is most 
evident at small radii, where magnetic diffusion has allowed the magne tic flux that would 



have been dragged onto the central object to occupy an extended region (ILi fc McKedll996 



Ciolek fc Konigll 1 19981 ; iTassis fc Mouschoviasl 120071 ) . The increased poloidal field strength 
tends to make the magnetic braking more efficient (jKrasnopolsky fc Konigll 120021 ) . Given 
these competing effects in opposite directions, it not at all obvious whether ambipolar dif- 
fusion would increase or decrease the overall efficiency of magnetic braking. Our detailed 
calculations have shown that, for realistic rates of cosmic ray ionization, the overall efficiency 
is increased, because the azimuthal field slippage is small and the poloidal field trapped at 
small radii is strong. The efficiency is reduced only for unrealistically low ionization rates, 
which speed up the azimuthal field slippage and lower the poloidal field strength by enabling 
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the trapped flux to escape to large distances. 



There is, however, ample evidence for large, rotationally supported disks on 10 2 AU 
scale or more around young stars, at least in Class I and II, and possibly in Class 0, sources. 
This begs the question: what is the origin of these disks? 

The most obvious possibility is that additional non-ideal MHD effects may weaken the 
mag netic braking fur th er. Ohmic dis s ipatio n is expected to be important at high densi- 
ties (IShu et al.l 120061 ) . iNakano et al.l (120021) estimates that it dominat es other processes 
at densities above ~ 10 12 cm -3 (see also iDesch fc Mouschoviad l200ll ). Since the den- 
sity in our standard model never exceeds this value, we believe that Ohmic dissipation 
would not significantly weaken the magnetic braking inside our computation domain (from 
10 14 cm< r < 2 x 10 17 cm). However, Ohmic dissipation may still enable a rotationally 
supported disk to form at a sm aller radius where the density is higher, as demonstrated ex- 
plicitly by lMachida et al.l (120081 ) using a resistive MHD core. They found a disk of ~ 10 12 cm 
in radius shortly (~ 11 days) after the formation of the so-called "second core"; the disk was 
not present in the ideal MHD counterpart. Such Ohmic dissipation-enabled disks can form 
inside our inner boundary (of 10 14 cm). Since the material crossing the inner boundary of 
our simulation is typically braked to a rotational speed well below the local Keplerian speed, 
any disk that may form later should be relatively small in the absence of angular momentum 
redistribution. 

Besides Ohmic dissipation, another non-ideal mechanism of magnetic flux diffusion is 
the Hall Effect. When dust grains abundances are low, the Hall Effect becomes important 
when collisions with neutrals decouple the ions (but not electrons) fro m the magnetic field. 
This occurs when the ion Hall parameter $ is less than unity (e.g., iNakand 1 19841 ). We 
have checked that > 1 everywhere inside our computational domain for the standard 
model. We do not expect the Hall Effect to change our results significantly, although this 
may be chan ged by the inclusion of dust grains. The exact effect depends on the grain size 
distribution (IWardle &: Ng!ll999l ). which is uncertain. Further study is required to determine 
the importance of the Hall Effect. If the additional non-ideal MHD effects do not allow the 
observed large scale disks to form, other means of weakening the magnetic braking must be 
sought. 



One possibility is protostellar outflow. Protostellar outflows are observed at all stages 
of st ar formation. They a re thought to clear away the envelope material as the source 
ages (lArce fc Sargentl 120061 ) . As the envelope is the anchor for the magnetic brake in our 
models, removing it may reduce the strength of braking sufficiently to allow large scale 
rotationally supported disks to form. This is particularly true for the more evolved class II 
and perhaps class I sources, where little envelope is left to brake the disk. It may not work 
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for class sources where t he outflows tend to be confined to the polar regions (e.g., HH 211, 



Gueth fc Guilloteaulll999l ). If this is the case, we would not expect these deeply embedded 
objects to harbor large scale, rotationally supported disks. They could still have relatively 
small disks enabled by, for example, Ohmic dissipation. It is difficult to determine the exact 
amount of rotational support in the equatorial region from current observations of Class 
sources. The situation will improve when ALMA comes online. 

We thank Mordicai Mac Low and Jeff Oishi for help in the implementation of ambipolar 
diffusion in ZEUS MP. This work is supported in part by NASA (NNG05GJ49G) and NSF 
(AST-0307368) grants. 
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Fig. 1. — The mass accretion rate of the standard model (solid line) compared to Shu's value 
for singular isothermal sphere (dashed). 
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Fig. 2. — Color map of density for the standard collapse solution, with velocity vectors 
showing simultaneous infall and outflow. The magnetic field lines (black contours) show 
strong pinching at small radii due to the collapse. The flattened high density region is a 
collapsing pseudodisk rather than a rotationally supported disk. 
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Fig. 3. — Radial (top panel) and rotational (bottom) velocities in units of the sound speed 
of the neutrals (thick) and ions (thin) on the equator for the standard model. Note the 
decrease of the rotation velocity with decreasing radius, indicating strong magnetic braking, 
especially in the deceleration region (shaded). 
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Fig. 4. — Velocity profiles of different collapse models. The top panel shows the neutral 
(solid) and ion (dashed) radial velocities in unit of the sound speed for three models with 
the same core mass-to-flux ratio A = 4 but different cosmic ray ionization rates £ = 10~ 16 
(blue), 10~ 17 (black), and 10~ 18 s" 1 (red). The rotation speeds of these models are shown in 
the second panel. The third and bottom panels show, respectively, the radial and rotation 
velocities of three models with a weaker magnetic field of A = 13.3 but different ionization 
rates C = 10~ 16 (blue), 1(T 17 (black), and KT 18 s" 1 (red). Note that the rotationally 
supported disk forms only in the model of the lowest ionization rate and weaker magnetic 
field. 



